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The dynamics of an elastic medium driven through a random medium by a small applied force 
is investigated in the low-temperature limit where quantum fluctuations dominate. The motion 
proceeds via tunneling of segments of the manifold through barriers whose size grows with decreasing 
driving force /. At zero temperature and in the limit of small drive, the average velocity has the form 
v oc exp[— const. /ft a / M ]. For strongly dissipative dynamics, there is a wide range of forces where 
the dissipation dominates and the velocity-force characteristics takes the form v oc exp[— S(f)/h], 
with S(f) oc 1/ f( d + 2 0/( 2 ~<) t ne action for a typical tunneling event, the force dependence being 
determined by the roughness exponent £ of the d- dimensional manifold. This result agrees with 
the one obtained via simple scaling considerations. Surprisingly, for asymptotically low forces or 
for the case when the massive dynamics is dominant, the resulting quantum creep law is not of 
the usual form with a rate proportional to exp[— S(f)/%]; rather we find v oc exp{-[S(f)/h} 2 } 
corresponding to a — 2 and fj, — 2(d + 2( — l)/(2 — £), with the naive scaling exponent for 
massive dynamics. Our analysis is based on the quasi-classical Langevin approximation with a noise 
obeying the quantum fluctuation-dissipation theorem. The many space and time scales involved 
in the dynamics are treated via a functional renormalization group analysis related to that used 
previously to treat the classical dynamics of such systems. Various potential difficulties with these 
approaches to the multi-scale dynamics — both classical and quantum — are raised and questions 
about the validity of the results are discussed. 

I. INTRODUCTION 

Static and dynamic properties of elastic systems in the presence of disorder have attracted the attention of physicists 
for more than three decades. Vortices in superconductorscl, charge density waves in solidsa, domain walls in magnetsa, 
and geological faultso are well-known examples of such systems. Mathematically, -some of the problems that arise in 
modeling are related to those in Burgers turbulenceO, stochastic growth of surfacesQ, and the stability of matterll The 
physics of dirty elastic systems thus impacts on several disciplines. In this paper we focus on the dynamic properties of 
driven elastic manifolds at temperatures low enough for quantum fluctuations to play an important role, in particular, 
the phenomenon of quantum creep. 

In many situations the elastic system can be driven by applied forces. For example, transport currents in supercon- 
ductors cause a Lorentz force to act on vortices, while a magnetic field applied to a ferromagnet produces an effective 
force on domain walls. One of the primary quantities of interest in such systems is the average velocity v of the 
manifold as a function of the applied force /, its dependence on the temperature, and on the magnitude of the ran- 
dom disorder that impedes the manifold's motion. The presence of the random pinning forces renders the theoretical 
analysis difficult, as perturbation theory breaks down in most of the important regimes due to the deformations of 
the manifold on many length- and time scales. However, some important general results are known: in the absence 
of thermal and quantum fluctuations the system stays pinned — i.e., the steady state velocity v — — if the driving 
force / is smaller than a certain critical value f c . At nonzero temperatures, even if the force / is smaller than f c , the 
system rwill move with a nonzero velocity due to thermal creep, see Fig. [l]. In the limit / <C f c the scaling theory of 
creepEnL£l predicts the law 

v cx exp [-U{f)/T] , (1) 

with U(f) — > oo as / — ► and there is no linear response to an applied force. The quantity U(f) can be interpreted 
as a barrier separating two neighboring minima of the free energy on some appropriate length scale that depends on 
the applied force; on long length- and time scales, the manifold's motion proceeds via thermally activated jumps of 
corresponding segments between different metastable configurations. 
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FIG. 1. Tunneling (lower dotted line) and activated motion (upper solid line) of an elastic string in a disordered potential 
under the action of a weak external force / <C f c - 



Several important assumptions are made in order to derive Eq. ([lj). First of all, at least in its simplest form, it 
is assumed that the barriers between different metastable configurations scale with length in the same way as the 
variations of the free energy that are caused by the balance .between the random pinning and elasticity; in general one 
could consider two different scaling laws for these quantities!!^. Second, Eq. ([!]) is an Arrhenius law with typical events 
dominating the creep-like motion. In general, rare events, in particular regions with anomalously high barriers, might 
lead to a modification of this law, perhaps to a form v oc exp [— (U(f)/T) ] with a ^ 1 an exponent characterizing the 
tails of the barrier probability distribution. Because of these assumptions, as well as the need for increased analytical 
understanding, it would be valuable to derive the classical creep law (|l|) starting from a microscopic description. 
Progress in this direction has recently been made, see Refs. (y],[l2[ where Eq. (|l]) has been derived from an at-least- 
partially controllable rcnormalization group (RG) expansion. Although these results support the validity of the scaling 
theory for thermal creep, serious questions remain which we will discuss later in this paper. We also note that neither 
the recent renormalization group calculations nor the present paper address an important aspect of creep in many 
contexts, in particular in the vortex lattice, i.e., the role of dislocations — the results are, so far, restricted to a truly 
elastic manifold -ssjith pinning weak enough that deformations are too small to induce dislocations or other strongly 
nonlinear effectsE^I. 

Experimental investigations have shown tJa&t even at very low temperatures creep-pf driven elastic systems, e.g., in 
vortex latticesEJ, of domain walls in magncts&3, and of charge-density waves in solidstj, still persists. This phenomenon 
suggests an explanation in terms of quantum tunneling-ni the manifold between different metastable configurations; 
the macroscopic manifestation of this is quantum creepEH. By analogy with thermal creep, one would guess that 

v<xwcp[-S(f)/h], (2) 

with S(f) a characteristic action for tunneling of a typical segment of the manifold whose size is determined by the 
applied force. In spite of a substantial number of studies of quantum creepta there is as yet no theoretical analysis 
that starts from a microscopic description of a driven elastic manifold interacting with impurities; the main goal of 
this paper is to present a first start at such an analysis, along with a discussion of the difficulties involved. We will 
show that a law of the form 

v<xexp{-[S(f)/n] a } (3) 

can be obtained using a renormalization group approach, as can the functional dependence of the action S on the 
external force / be found. Whether this is indeed the correct behavior is addressed at the end of the paper— 

A natural tool to investigate quantum creep is the functional renormalization group (FRG) expansionlla, an e- 
expansion near the upper critical dimension (d c = 4) of the random manifold problem. This involves, intrinsically, 
an infinite number of marginal operators that can be combined into one or more functions. We will find that in this 
approach a natural separation of the frequency scales occurs between inter-valley and intra-valley fluctuations of the 
manifold. The inter-valley fluctuations correspond to the motion of the manifold on large scales on which there are 
many separate valleys in the "landscape" caused by the randomness. The intra-valley fluctuations correspond to the 
much faster motion of the manifold within one valley. 

The paper is organized as follows. In section III we formulate the model and find the appropriate effective action. 



In section III we show how to derive the creep law (|2|) using scaling arguments. Next, we study the problem using 
the renormalization group expansion (section |Tv| ) and then summarize and discuss the main results. We end with 
an analysis of the limitations and problematic aspects of both our and previous approaches to the creep problem in 
section [V| 
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II. MODEL AND EFFECTIVE ACTION 



We describe the elastic medium by a non-dispersive elastic manifold interacting with a quenched random potential; 
this generic model captures all the essential physics while allowing for an extensive analytical treatment. The classical 
motion of the displacements u(z, t) of the manifold away from its undeformed state is described by the dynamical 
equation 

V d t u + pd 2 t u = cW 2 u + F{u, z) + / + f th (z, i), (4) 

where the friction force rjdtu and the inertia pd^u are balanced by the elastic- (cAu), random pinning- (F(u,z)), 
driving- (/), and thermal- (/ t h(z,i)) forces. We will consider the case, applicable to domain walls, of d-dimensional 
manifolds in d + 1-dimensional random environments, so that z £ TZ d and u £ 1Z 1 . 

The random pinning force F(u, z) is taken to be Gaussian random with mean zero and a correlator 



F(u, z)F{u', z') = A(u - u') 5{z - z'), (5) 

where the overbar denotes an average over the disorder. The function A(u) decays rapidly with a characteristic scale 
£. The thermal noise / t h(z,i) is Gaussian white with a correlator 

</th(z, £)/th(z', 0) = 2vT6(t - t')5(z - z') = K (t - t') 5(z - z'), (6) 

with (. . .) denoting the average over thermal fluctuations. 

In classical statistieaJ-^nechanics we can reformulate the problem posed by a stochastic differential equation as an 
effective field theorytSEJ with the help of the Martin-Siggia-Rose (MSR) approach. After averaging over thermal 
fluctuations the MSR-action corresponding to (Q) has the form 

^4msr = — J d d zdt iy{z, t) (jidtu + pd^u — cAu) 
d d zdtiy(f + F(u,z)) 

+ \ j ddzdtdt ' iyi 2 ^)^ ~ 0%( z >*')) ( 7 ) 

with y(z, t) an auxiliary field used to enforce the equation of motion. The probability of a particular dynamical 
evolution w(z,i) under the stochastic process is proportional to J T>[y] exp{AMSR[w, y]}. 

The main goal of this paper is to investigate the influence of quantum fluctuations on the system whose classical 
limit is described by Eq. (|4|). Again, it is convenient to formulate the problem as a field-theory, i.e., to write the 
effective action describing the elastic system in the presence of quantum fluctuations in a way analogous to (|). When 
calculating time-independent quantities in equilibrium systems (i.e., in the case / = 0) the corresponding quantum 
action has a Euclidean form and the quantum partition function can be written as an imaginary time path integralEil. 
However, here we are interested in non-equilibrium properties (/ ^ 0) and the Euclidean action cannot be used. 

A fopaalism allowing to study the real time dissipative quantum mechanics of a system is that of Feynman and 
VernonE3: the quantum amplitude ^(xf,tf,Xi,ti) for a system to have a coordinate Xf at time tf, if at time ti it had a 
coordinate Xi, can be written as a path integral f Xf l f D[x\ exp(iS[x}/%), with S[x] the classical action. Consequently, 
the probability of the transition Xi,U — > Xf,tf can be written in the form 

rx f ,tf px f ,t f 

V{x f ,t f ;x t ,ti)= I I V{x}V[x'}exp(iS[x]/h)exp(-iS[x'}/h). (8) 

J Xi ,ti J Xi ,ti 

In non-equilibrium quantum mechanics the functional @) plays the same role as the partition function in equilibrium 
problems. We see that the effective action can be written as iS[x] — iS[x') and includes two different paths [x(t)] and 
[x'(t% 

The standard way to include dissipation is to write the action S[x) in the form S[x] — Sq[u] + S'bath [-^bath] + S'intM 
representing the actions corresponding to the elastic manifold, the bath, and the interaction between the bath and the 
manifold, respectively (note that x = (u, A"b a th|)-and x' = (u', Aj^ath) are the coordinates describing both the manifold 
and the bath). In the Caldeira-Leggett modeltJ the terms S'bath [^"bath] and 5*i n t[a;] are taken to be quadratic in the 
bath coordinate Abath and hence can be integrated out at the expense of effective interactions that couple different 
times. The resulting action takes the form (see appendix [A] or Ref. |24|) 
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iS[u, y] 



~ J d d zdt y(z, t) (j]d t u + pd 2 u — cAi 
d d zdtdt' iy(z, t)n(t - t')iy(z, t') 



1 



- I d d zdl 



(9) 



with 



and the random potential 



u(z,t) = (u(z,t)+u'(z,t))/2, 
y(z,t) = u(z,t) - u'(z,t), 



U(u,z) = - 

The Fourier transform of the effective noise correlator n(t- 
dissipation theorem, 



r 

I duiF{u\,z). 



■t') of the bath fluctuations obeys the quantum fluctuation- 



k(uj) = rjhu> coth 



ha; 
2T' 



(10) 



where r\ is related to the spectral density of the bath, see appendix [A|. 

The action (^) allows one, in principle, to investigate all the properties of a driven elastic medium in the presence of 
quantum fluctuations. The classical limit, the MSR action ^4msr = iS/Ti can be recovered by expanding the random 
potential energy U in tiie y-field and substituting u by u and y/Ti by y. Such an expansion is valid at high temperatures 
and it has been argucdn3 that it also applies to the limit h — ► at T = 0. This suggests that the dynamics of quantum 
systems in the semiclassical regime, where tunneling processes through large barriers dominate, can be described by 
Eq. ([!]) adth the noise obeying the fluctuation-dissipation theorem (|l0|); this is called the quasi-classical Langevin 
equationBj (QLE) approach. In a number of paperso'LJ the QLE has been used to study the non-equilibrium dynamics 
of quantum systems. In particular, in Ref. ^6| the quantum tunneling of an overdamped quantum particle has been 
studied using this approximation. Writing the inverseriifetime T of a metastable state in the form T = Pe~ s / h with S 
the tunneling action and P the prefactor, it turns outcJ, that the QLE gives the correct value of the tunneling action 
S up to a multiplicative factor of order unity (note that for quadratic systems the QLE is actually exact, as is the 
expansion U(u + y/2, z) — U(u — y/2, z) w U'(u, z) y for this case). 

In this paper we use the QLE to study quantum creep in a disordered medium. Although we would like to do 
better, we are primarily interested in obtaining the correct scaling laws and asymptotic forms; the exact coefficients 
are of less interest. The action ^ with the quantum noise ( |l0| ) is much simpler to analyze than the full quantum 
action (^|). We believe that the QLE should be a useful tool for investigating the dynamics of quantum systems 
as it appears to capture the essential physics, while being — at least relatively — tractable analytically. Quantum 
fluctuations appear in the QLE approach as an effective random force with a correlator n(t), see (^0|), acting on the 
system. The quantum |decay of a metastable state then is equivalent to the thermally activated escape of a particle 
driven by colored noisec3. 



III. SCALING ARGUMENTS 



A. Statics and Thermal Creep 



An elastic manifold subject to a random potential behaves as an assembly of approximately independently pinned 
segments of a characteristic size that is determined by the balance between the pinning and elastic forces. The driven 
motion of the manifold is dominated, both classically and quantum mechanically, by successive jumps of segments 
of the manifold between subsequent metastable configurations. In order to understand this behavior one needs to 
consider various characteristic scales. n 

The characteristic size of the collectively pinned segments is of ordertl L c ~ [c 2 ^ 2 /A(0)] 1 ^ 4 ~ d - ) , the Larkin length 
over which the typical distortion of the manifold \u(L c ) — u(0)\ is of order £, the scale of the correlations in the 
random pinning landscape. The characteristic energy scale of the deformations on scale L c is U c ~ c^ 2 L d T 2 . Beyond 
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L c , the typical displacements of the manifold grow as \u(L) — u(0)| ~ £,(L/L C ) ( ' with £ the wandering exponent, 
£ = (4— d)/3 for the random force case, while (d,i ~ 0.2083 (4— d) for a short-range correlated random potential; the 
difference in energy between metastable configurations deviating on a scale L by a displacement Su(L) is controlled 
by the balance between the elastic- and pinning energies, both of which have a magnitude of order U C (L/L C ) with 
the energy exponent 6 = d + 2Q — 2 and the wandering exponent C given above. 

In the presence of an applied force / it is favorable for a segment of size L to move into the neighboring metastable 
valley if the energy fL d Su(L) ~ / 'L d ^(L/ 'L c ) gained due to the presence of the force is larger than the difference 
of the elastic- plus pinning energies between the two configurations; for small forces, this will only occur for large 
segments. Equating the two energy scales, we obtain the minimum characteristic size of a segment that can move to 
the lower energy sitctj, 

/ f \ 1/(2-0 



/ 

where f c ~ c£/L;? is the critical applied force needed to move the manifold in the absence of thermal or quantum 
fluctuations. The scale Lf is the minimum size of a segment that can hop in order to lower its energy; the macroscopic 
motion then proceeds via jumps of segments of this characteristic size over distances 

/ f \ C/(2-C) 

»'~<7) (12) 

separating two neighboring metastable configurations. 

Consider first the classical motion, i.e., thermal creep: The conventional assumption is that the height of the energy 
barrier that must be surmounted for such a motion to occur classically is determined by the scaling of the static 
energy U C (L/L C ) 6 (we will discuss the validity of this assumption later in this section). The characteristic height of 
the barrier that dominates the motion at a small force / then is 

U(f) ~ cu}L d f 2 ~ U c (L f /L c ) d+2 <- 2 ~ U c (fc/f) {d+2C - 2)/i2 ' C) ; (13) 
this yields an average velocity for thermal creep, see Ref. 

v ~ exp[-t/(/)/T]. (14) 



B. Quantum Creep with Dissipative Dynamics 

The quantum tunneling of segments of the manifold can be estimated in a manner analogous to their classical 
thermal activation: Let us assume that the manifold is located in one of the metastable configurations and subject 
to a small driving force /. One expects the lifetime of such a state to be proportional to exp(S/h), with S the 
characteristic action describing the tunneling through the barrier separating the metastable minimum from one with 
a lower energy. We estimate S using the standard theory of the decay of metastable states. We first assume that the 
effects of dissipation dominate over those of inertia; the Euclidean action of the manifold then can be written in the 
form 



+°° _ . ., +oo . _ „ 



'EuclN 



dT / d d 



du Y i jj / t , V f ,(u(z,t)-u(z,t')) 



2-UJ +U ^~f U+ t J ^ (r - r0 



>2 



(15) 



with r the imaginary time and rj the dissipative coefficient. In order to find S one needs to find the imaginary time 
trajectory connecting the initial and final configurations of the manifold and calculate its action. We emphasize that 
as our problem is a non-equilibrium one, the full quantum action (Q) should be used. However, at low driving forces, 
the problem we study is a quasi- stationary one, as the lifetime of a metastable state is large and the manifold can 
be considered to be in local equilibrium. Therefore we can use the simpler (Euclidean) action (|l5| ) instead of the full 
dynamic action (0). 
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/ U(f) 
















FIG. 2. Tunneling of a particle with an effective mass pLf and dissipative coefficient r}Lf through a barrier of height U(f) 
and width Lf. 



In estimating the tunneling action we find the size Lf of the tunneling segment as in the classical case of thermal 
creep, since this is determined by the static balance between the variations in the pinning energy and the energy gain 
due to the external force. Furthermore, we assume that there is only one characteristic time scale rj associated with 
the tunneling process. We can then consider the tunneling event, crudely, as an effective process in which a point-like 
object with a friction coefficient r\f = rjL'j tunnels through a potential barrier of height U (/) and extent u / , see Fig. |^. 
The saddle-point of the action ( |l5| ) is minimized when 



U(f) 



Vf — 

T f u f 



(16) 



where we have substituted the dynamic term by its characteristic value and U(f)/uf is the characteristic force acting 
on the "particle". The tunneling time 77 is given by r/ ~ rjfU 2 f/U(f), resulting in the tunneling action 



S(f) ~ s ri 



j,(d+2Q/(2-C) 
7) 



(17) 



with S v ~ il£, 2 Li the characteristic scale of the dissipative action of a segment with the diameter of the order of the 
Larkin length L c . The velocity of the manifold is related to S(f) via (||). 

The characteristic crossover temperature T cr from quantum to classical creep can be obtained by comparing U (/) /T 
and S(f)/U, i.e., 



T 



nu{f) hu c 



2/(2-0 



(18) 



jives the scaling conjecture for the dissipative limit. 



C. Quantum Creep with Inertial Dynamics 

We now assume that the inertial term in the action is more important for the tunneling process than the dissipation. 
In this case the Euclidean action analogous to that given by (Il5|) can be written in the form 



+00 



SeuciM 



dr 



d d z 



du 



+ U(u, z) - fu + 



du 
l)t 



(19) 



with p the mass density of the manifold. The size Lf of the tunneling segment is the same as in the dissipative case. 
In order to find the action S( f) one again needs to find the tunneling time 77; comparing kinetic and pinning energies 
for the saddle-point configuration, we find (with pf 



Uf 



pLj the effective mass of the tunneling object) 
U(f) 



(20) 



G 



1/2 

i.e., Tf ~ [pf/U(f)) Uf, and the tunneling action can be estimated as 

(d+2C-l)/(2-C) 

with S p ~ (pUcL^) 1 / 2 ^ the chara cteris tic action with inertial dynamics on the Larkin scale L c . The quantities L c and 



S(f)~S p (fy) , (21) 



U c have been defined in section HI A above. The characteristic crossover temperature can be determined as in the 
dissipative case and is given by T cr ~ (HU C / S p )(f / f c ) • 

D. Assumptions 

Two important assumptions have been made in order to derive the scaling prediction (ph: First of all, it has been 
assumed that the barriers for the manifold's motion scale in the same way as the variations of the energy. If there are 
two independent exponents 9 = d — 2 + 2£ parametrizing the static energy scaling and ip > 9 characterizing the scaling 
of energy barriersr 9 H 30 l, then the barrier for thermal creep will have the form U(f) oc L^. The simplest assumption for 
the quantum motion is that this same barrier height is the appropriate one for the quantum barrier traversal. While 
the dissipative result ( p7| ) would remain the same (as S^iss ~ Vf u2 f does not depend on the height of the tunneling 

barrier), the result for the inertial dynamics would take the form 5i ne rtiai oc Ly+^+zC)/ 2 _ 

We should emphasize, however, that if ip > 9 then this is because the dynamics on intermediate length scales affects 
that on the large scale Lf. A dependence of the large scale dynamics, of interest for quantum creep at small applied 
forces, on the intermediate scale physics can also occur quantum mechanically; but the way in which it occurs could be 
quite different from the classical situation and it is possible that each type of quantum dynamics has its own exponent 
tf)Q distinct from the classical value ijxA. Once we admit the possibility of a dependence of the tunneling dynamics of 
large segments on the dynamics at intermediate scales, the basic assumptions of the scaling arguments given above 
break down. But this is exactly the kind of problem for which a renormalization group framework is needed. 

Another potentially important effect that was not taken into account in the simple scaling analysis summarized 
above are rare events: in general, it is possible that the velocity is not controlled by the typical barriers — or, 
more precisely, by the tunneling through these typical barriers — but by rare anomalously large barriers. In the 
simplest scenario, this could give rise to a quantum creep law of the form v oc exp{— [S(f)/h] }, with a =/= 1 a 
non-trivial exponent characterizing the tail of the distribution of barrier heights. Just such a phenomenon controls 
the temperature dependence of the diffusion coefficient of a single particle diffusing in a Gaussian random potential^, 
see the discussion in section [V E| . One mechanism for an exponent a > 1 is that rare, rather than typical, events 



dominate the macroscopic motion. 

IV. RENORMALIZATION GROUP ANALYSIS 

A. Derivation of RG Equations 

Having presented the basic scaling arguments for both classical and quantum creep, we now turn to the main task 
of this paper, an attempt to derive the creep laws from a microscopic description of the dynamics. With the roughness 
exponent £ vanishing in dimensions d > 4, one expects that d c = 4 is the upper critical dimension of randomly pinned 
clastic manifolds, an expectation which turns out to be correct both for the equilibrium properties in the absence of 
a driving force and for the fluctuationless (T = 0) critical depinning transition at f c in the presence of a driving force 
/. A renormalization group e-expansion about d c = 4 then allows one to investigate the various large scale properties 
of the pinned elastic manifoldES. An important feature of this renormalization group analysis is the fact that even 
to lowest order one has to take into account infinitely many variables, in practice by renormalizing a function that is 



essentially the correlator A[u(z) — u(z')] of the random force (but see the discussion in Ref. |33| and in section VB ) 

The functional renormalization group (FRG). approach has been successfully applied to both static equilibriur 
and fluctuatioaless (T — 0) depinning problemsOEi Recently, the FRG has also been applied to the thermal creep 
problemy'tSEa, although its applicability is far more problematic. For forces / <C f c , the motion proceeds either by 
activation or tunneling of large segments between different metastable minima and thus the dynamics is intrinsically 
non-perturbative: in particular, the dynamic exponent z, which is shifted only slightly from its naive value of two near 
the critical depinning transition, must be radically modified in the creep regime. Nevertheless, it appears (though, see 
section [V]) that for thermal depinning one can handle this within the FRG and renormalize until a scale is reached 
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beyond which the disorder can be taken into account perturbatively and the equation of motion can be solved. By 
retracing back the FRG flow, one can then determine the average velocity of the manifold. We will follow this approach 
here, but return later in the paper to question its validity both for the quantum and classical cases. 

Formally, in both classical and quantum cases, the RG analysis of creep involves the flow away from the static 
(/ = 0) fixed-point that controls the undriven pinned system under the action of a small applied force /; under the 
RG transformation the parameters describing the system flow in such a way that eventually the effects of disorder can 
be neglected beyond a certain length scale. But there is an important difference between the classical and quantum 
cases: In the classical casalj one can restrict the analysis to the low-frequency limit, i.e., including the friction term 
(rjdtu) alone into the renormalization procedure is sufficient. This is a consequence of the fact that, to exponential 
accuracy, classical creep does not depend on the dynamics (e.g., inertial versus dissipative) of the system. In the 
quantum case, it will turn out that an analogous treatment leads to a spurious "loca lizatio n transition" where the 
average velocity of the manifold drops to zero at a finite length scale ~ L c , see section [IV B . In order to carry out a 



correct analysis one must include all the frequencies into the renormalization group; specifically, we must make the 
substitution 

/t °° 
dTD{r)u{t -t) = r,d t u + ^ rt {k) d*u. (22) 

-°° k=2 

It will turn out that, even if initially rp k ^ = for all k > 2, the dynamic parameters rp k ^ will grow rapidly under the 
flow and become crucial; the renormalized spectrum Di(u) is very different from the bare linear spectrum —iu>r). 

It has been argued above that in order to investigate quantum creep we can use the action (^) with the correlator 
given by jlO|). It is convenient to perform the disorder average in (Q) and work in a frame moving at the average 
velocity of the manifold, i.e., we substitute u(z,t) — > u(z, t) + vt. We also introduce the above general dynamics D; 
the resulting action, averaged over both quantum fluctuations and random pinning, takes the form 

A = - J d d zdtiy{t)®u- cd^u) + i J d d zdtdt' ' iy(z,t)n(t - t')iy(z,t') (23) 

+ J d d zdtiy(f - rjv) + ^ / d d zdtdt' iy(z, t) A[u(z, t) - u(z, t') + v(t - t')] iy(z,t'). 

The correlation functions corresponding to the quadratic action in the absence of pinning are given by 

u(k,uj)u(k',uj') = ^ J( y(w + </)5(k + k') = C(k, lo)S(uj + u/)6(k + k'), (24) 

|ck 2 + D(u>)\ 

1 

ck 2 +~D(lu)~ 



u(k, u)iy(kf, oj') = —5 =rr^S(u + u')6Ck + k') = R(k, oj)5(uj + o/Wk + k'), (25) 



with the usual definition of the Fourier transform 



/(k,w) = J d d zdtf(z,t)cxp(-ikz + iu;t) (26) 
and D(lu) sa r)(—iuj) for small u>. Because of the fluctuation dissipation theorem, we have 

k(cj) = -MmD(u) coth — . (27) 

In the regime of interest for quantum creep, the motion of the elastic manifold can be represented as jumps between 
different metastable states and the lifetime of each metastable state is large. It is thus possible to extract some 
information about the dynamic properties of the system, in particular, the average velocity v, from static equations 
if we cutoff the renormalization group flow at the relevant length scale, here, Lf ~ L c (f c / /) 1 ^ 2 ~^', see section [II (a 



similar situation arises in the conventional theory of the decay of metastable states: although the original problem is 
formally non-equilibrium, the system is trapped for a long time in the local equilibrium state corresponding to a local 
minimum of the free energy; one can approximate the partition function of the system by that calculated assuming 
quasi-equilibrium and then find the flow of the probability out of the local potential well). This is why we could use 



the Euclidean action instead of the full quantum mechanical action in section III. For values of the RG variable I 
smaller than ln(L t / L c ) we then can use the appropriate static versions of RG equations. 

The RG flow equations are obtained by integrating over fast modes in the effective action ( p3| ) with the last two 
terms in the action treated as perturbations about the quadratic action. We write all the equations, except those for 
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the correlator A; (it), the temperature XJ, and the Planck constant hi, at finite velocity. We will use the subscript 
I in order to distinguish renormalized quantities at length scale L = A _1 e ( from bare quantities without subscripts, 
e.g., r/ = r/o, T = To, h = Hq, A(it) = Ao(it), etc. Introducing the large momentum cutoff A we obtain to leading 
non-trivial order in A;, 

diA t (u) = (4-d-2()Ai(u) + (uA' l (u) + C>{A,t = 0)A'{(u) 

+ IA'{{u) [A z (0) - Aj(u)] - IA; 2 ( m ), (28) 

dji = (2 - 2C)/, + / dtRf^t^vtt), (29) 



ft«l = («-0»l, (30) 
= (4-d-2C)/si(t) + ;rtft/ej(t)-C^(A,t)A{>jt), (31) 

ftA(w) = 2D,M " zwdMw) - [ ^-A'{ (uj'/vi) [R?(A,u>') - R>(A,u +cj')] , (32) 

J 27T Vl 

diTi = (2-d-2QT h (33) 
dihi = (2-d-2t-z)h h (34) 

with / = A(iA d /c 2 A 4 defined in terms of the surface area Ad of a d-dimensional unit sphere. Both the dynamic 
exponent z = z(l) and the roughness exponent £(Z) are at our disposal to adjust for convenience; it will generally be 
most useful to choose £ to be the /-independent value giving rise to a well behaved fixed point function A* (it) in the 
absence of fluctuations or drive. How best to adjust z(l) we reserve for later; conventionally it would be chosen to fix 
the coefficient in the low- frequency part of Di(u), in the dissipative situation studied here, of rji. The shell- restricted 
correlation functions C > (A,t) and i? > (A,t) are given by 

c ' >(A ' tH t$/^"' c < A '"» < 35 » 

and 

R>(A,t) = ^ [ pe-^R(A^); (36) 

here, C(A,o>) and R(A,u>) are defined via ( p4| ) and ( p5[ ) with D(u>) substituted by Di(m). Also, we define the coefficient 

r, = C>(A,i = 0), (37) 

which will appear frequently in the following. We chose an initial dynamic spectrum of the form 

Dq(lj) = —ii]uj + pu 2 . (38) 

The excess of the applied force over that needed to sustain the motion of an unpinned system with dissipative 
coefficient rji is 

.fi = fi-Vivi- (39) 

Above, we have defined the Fourier transform of the second derivative A" (it) = d 2 A by A"(p) = J due lpu A"(u), so 
that 

/+oo 
dte lut A'({vt)- (40) 
-oo 

this quantity plays a crucial role in the dynamic reeprmalization. Eqs. (|28]), (p9|), (^l|), (^0|), and ( |33| ) have been 
obtained before in the discussion of classical creep!3; they are the same in the quantum case, while Eq. d32| ) is 
different, however. We now show how it can be obtained and in what respect it differs from its classical analog. 

After averaging the disorder term over the fast modes one obtains the following feedback SA to the term 
— J d d zdtiy(z, t)Di ® u in the action, 

5A = - [ d D zdtdt' iy(z, t)R> (A, t - t') [u(z, t) - u(z, t')] A'{ [v(t - t')} dl. (41) 
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Writing the functions A; far) and i? / > (A,r) as Fourier integrals we obtain Eq. (p2[), which explicitly includes the full 
dependence of the displacement field u(z, t) on time, i.e., it is nonperturbative in frequency. In the classical analysis of 
Rcf. [[2] it was assumed that the low- frequency limit of the RG equations is valid for all frequencies, i.e., for u — > the 
only important term is that containing r/idtu (corresponding to —rjiiuju{Lo) in Fourier space). The equation analogous 
to ( |32| ) then is given by its low-frequency limit, i.e., instead of the full dependence u(z, i) — u(z, t') one considers only 
the first term of its Taylor expansion in t — t'. It was also assumed in Ref. [l^ that the response function -R^A, t) and 
the correlation function C(" (A, t) depend only on the frictional coefficient r\i . In the classical limit this approach leads 
to reasonable results. By contrast, ignoring the frequency dependence of the dynamics in the quantum case leads to a 
spurious "localization transition" , implying a zero average velocity of the manifold below a small but n on-zero driving 



force / and in the presence of quantum fluctuations K ^ 0. We will discuss this issue in section IV B 



In the case of a purely dissipative dynamics there is another problem that occurs even at the initial stage of the 
renormalization: the integral in the expression for C^(A,t — 0) diverges at large frequencies, see Eqs. ( |35| ) and (|24|). 
For large u> the noise ki(u>) (see (jlG)) and the dynamic spectrum Di{uj) are both proportional to to for a dissipative 
dynamics and, consequently, Cf'lT^t = 0) diverges logarithmically as J dto/uj at large frequencies, see (pij). This is 
not unexpected: it is just such a logarithmic frequency dependencej-that can cause localization in models with a single 
degree of freedom coupled to a bath that provides a linear frictiorO; thus integrating out all the frequencies at once 
is problematic even for short wavelength deformations of the manifold. In reality, we expect that an inertial term 
pdfu or some other frequency dependence describing the small scale dynamics will provide a cutoff at high frequencies 
(at Wo ~ f]/p f° r the inertial case); alternatively one could introduce a sharp cutoff by hand. We will consider both 
possibilities later, noting now that how this is done will affect the results far more than one might expect. 



B. Structure of the Quantum RG Flow 



The main goal of the remainder of this section is to analyze the system of RG equations derived in section IV A in 
the limit of small driving forces. We first discuss the important features and then provide a more detailed analysis 
in the following section. First, let us analyze the renormalization ( p8| ) of the force- force correlator. In the absence of 
quantum and thermal fluctuations (i.e., k(ui) = for all u>), the correlation function C ; > (A,t = 0) is zero. In this case, 
the nonlinearities in the flow equation of the function A; fa) cause it to become nonanalytic on a finite length scale, 
even if the bare correlator A fa) is analytic, see Ref. [l8[ This is easily seen by differentiating Eq. (E8|) twice with 
respect to u and substituting u = 0, resulting in an equation for the evolution of the quantity A"(0) alone. The simple 
autonomous flow equation for A"(0) leads to a divergence at a finite scale l c ~ e" 1 ln[c 2 /A d ~ 4 A '(0)], producing the 
Larkin- or pinning length L c ~ [c 2 /A"(0)] 1 ^ 4_d - ) ~ [c 2 ^ 2 /A(0)] 1 ^ 4_d ^ ) where collective pinning goes over into strong 
pinning. Beyond L Cl the perturbative description breaks down as multiple competing minima appear in the pinning 
energy landscape. The infinite second derivative suggests that the function A; fa) will have a discontinuous first 
derivative at u = 0, i.e., AJ(+0) = — AJ(— 0) 7^ 0. On length scales shorter that L c , the smoothness of A/ fa) reflects 
the smooth reversible evolution of a segment as it is pulled by its neighboring regions. But on larger length scales 
the internal deformations of the segment will cause it to jump discontinuously and irreversibly from one metastable 
configuration to another — this is reflected in the discontinuity in AJ . At the scale L c the force correlator essentially 
has reached its fixed point shape with a height A*(0) ~ (cA 2 ) 2 £ 2 A -rf e~ 2 ^ c and a width £* ~ £e~^' c , resulting in a 
cusp |A*'(0+)| - (cA 2 ) 2 £A- d e - c '<=, see appendix 

ic.tuations the coefficient Ti = Cf~ (A, t = 0) becomes nonzero and 
The derivative AJ(0) of the correlator at the origin is zero for 



In the presence of small quantum or thermal J. 



this leads to a smearing of the cusp in A;fa)t! 
C > (A,t = 0) 7^ but changes rapidly in a boundary layer around the origin. The cusp_that was present in the 
absence of fluctuations is smeared over a region ii sm oar which can be estimated by comparing^ the terms Ti | A"(0)| ~ 
^ ^smear/j /^smear and [A d /(cA 2 ) 2 ]Affa) in @. The derivative AJfa > it sm0 ar) approaches the fixed point 
value A*'(+0) found in the absence of fluctuations and we obtain the boundary width it smoa r ~ ri(cA 2 ) 2 /A d |A*'(0+)|. 
The curvature A"(0) ~ — [A d /(cA 2 ) 2 ] | A*'(0+)| 2 /Ti then diverges as the fluctuations renormalize to zero on large 
scales. 

Physically, u S mear(0 can be understood in terms of the equilibrium response of a segment of size L = K^e 1 to 
the motion of its neighboring regions. Usually, a small displacement of neighboring regions will cause only a small 
readjustment of the segment of interest within its local energy minimum. The exceptions to this occur when the 
minimum energy configuration of the segment jumps from one configuration to another as its neighboring regions are 
slightly displaced: it is these jumps that give rise to the cusp in A fa) in the absence of fluctuations. At non-zero 
temperature or in the presence of quantum fluctuations, the behavior will not change much except near these jumps, 
where there will be a range of positions of the neighboring regions over which the segment of concern will have a 
non-negligible probability to be in either of two distinct configurations. This will result in a smearing of the cusp 
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in A; (it) over the scale u SIaeax (l) of neighboring region displacements on which this split probability typically occurs. 
On large scales, the fact that the energy scale grows as L e means that it is much less likely that the position of 
a segment will fluctuate between two energy minima. This is reflected in the renormalization towards zero of F; 
and the concomitant flow towards zero of the smearing scale it sm ear ^ F;. Summarizing, the function A; (it), whose 



renormalization is given through equation (28), develops under the RG in the following way: At the scale L c , A; (it) 
is close to its fluctuationless fixed-point function A* (it), with fluctuations affecting its behavior only in the vicinity of 
the point u — via a smearing of the discontinuity in A[ in a boundary layer whose size is controlled by Ti . Once we 
know how the function A; (it) evolves, we can substitute it into the other RG equations and see how other quantities 
renormalize under the RG flow. 

In addition to u smcar , there are two other important displacement scales: the characteristic scale £* of the fixed- 
point correlator A*(it) and the scale u vc \ associated with the velocity it vo i ~ ?y/V//cA 2 . The latter is the product of 
the velocity vi and the characteristic timescale rji/ch? of the low-frequency part of the response function i? / > (A, t) at 
wavelengths of the order of the cutoff A -1 . At scales somewhat larger than the Larkin length L c cx A _1 e' c , 

u ve l -C u smeax -C £*. (42) 

This is a consequence of the fact that vi is exponentially small (cx exp{— [S(f)/h] a }) in h, while itsmear proportional 
to a power of hi and £* is a static quantity that does not depend on hi (in the classical case the relevant displacement 
scales obey the same relation with the role of h played by the temperature T). During the RG flow, u smG ar decreases 
gradually with decreasing hi whereas u vo i increases rapidly due to the sharp increase in the viscosity i]i . 

Eventually at some scale both u smG ar and it vc i are of the same order. At this scale, which is the crossover scale 
Lf appearing in the scaling arguments for creep, the functions A[(vt) and A"(vt) are roughly given by A* (+0) and 
A*"(+0) on the time scales that dominate in Eqs. (p9[), (|3l|), and (p2|). At even longer scales, u ve i eventually becomes 
comparable to £* , see Rcf . [l2|. Beyond this scale the randomness is effectively smoothed by the motion of the manifold 
and can be neglected or treated perturbatively. 

In analyzing the renormalization of other quantities by the random pinning forces, we thus see that there are two 
important regimes: 

i) Uvci < M smcar , the nucleation regime, (43) 

where we can approximate AUvt) by AJ(0) = and A"(vt) by A"(0). And 

ii) "smear < w ve i < £*, the depinning regime; (44) 

it was arguedEl that this latter regime resembles that of the critical depinning transition!! in the absence of fluctu- 
ations; we can then approximate AJ(ui) by A*'(+0) and A"(vt) by A*"(+0). These two regimes are separated by 



the scale Lf discussed in section HI, We will now show how the scale Lf ~ £c(/c//) 2 ® naturally appears in the 
solution of the RG equations and separates the two regimes. 

In order to do so we concentrate on the force equation (|2^) which involves the slope A[(vit) of the disorder correlator 



near the origin. The latter is small and vanishes at u = for I < l c , hence // = e^ 2 grows exponentially (since the 
velocity vi of the manifold is exponentially small we can neglect the term rjivi in the equation for //). In the absence 
of fluctuations the slope A[(vit) ~ A*'(0+) rapidly turns on as the correlator forms a cusp at l c ; if the disorder term 
(A d /cA 2 )|A*'(0+)| overcompensates the scaling term (2 — C,)e^ 2 ~^ lc f, the force will start renormalizing to zero while 
in the opposite case it will continue to increase. Balancing the two terms we find the critical force density f c ~ c£/L 2 . 

Fluctuations smearing the correlator on the scale « snlcar soften and delay the depinning transition until the growing 
width it ve i of the response function encloses the emerging cusp in Ai(u) at If, u ve \(lf) ~ u S mear('/)- Below If we 
have A[(vit) ~ A[(0) = and the force increases exponentially // = e^ 2- ^/; starting with a small force / -C / c , 
depinning occurs as the cusp emerges at u smcal (lf) ~ u vc \(lf). Replacing again AJ(v;t) ps A*'(0+) in (|29| ) we obtain 
the depinning condition //, ~ (A d /cA 2 )|A*'(0+)| ~ / c e^ 2_ ^-" c ; the depinning scale If then relates to the force / via 
(2 - 0(1/ - U ~ Hfc/f) or 

L f ~L c (fc/f) 1/(2 - C) . (45) 

We thus see that the characteristic length scale for creep — the size Lf of the segments that can jump to lower the 
energy — appears naturally within the RG framework. We emphasize, however, that Lf is entirely determined by 
static properties; so far we have not had to address the crucial issue of dynamic renormalizations — to these we now 
turn. 

The crucial quantity needed for the dynamic renormalization is A"(0), which can be, obtained from Eq. ( p8[ ) by 
comparing r;A"(u) with the nonlinear terms evaluated at the origin it = 0. This yieldsl3 
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(46) 



(see appendix |b|) and enables us to find the renormalization of the viscosity r\\ . Since we are primarily interested in 
calculating the velocity to exponential accuracy, the calculations simplify significantly. If we do not renormalize time 
scales explicitly, i.e., choosing z(l) — 0, the velocity v is given approximately by the viscosity at the scale Lf, 

In In — . (47) 

The validity of this condition to exponential accuracy is due to the following argument: At the scale Ldyn where 
Uvd ~ r]ivi/cA 2 ~ £* we have a velocity vi cx l/r]i as the effects of pinning are negligible beyond this scale. On 
intermediate scales, i.e., between Lf and Tdyn both rji and vi renormalize but not exponentially; the only exponentially 
large renormalization originates from scales between L c and Lf. Thus, to exponential accuracy, we can ignore the 
renormalization on scales larger than Lf (and smaller than L c ) and justify ( p7| ) (note that, in general, only such 
dynamic quantities as rji and vi are exponentially renormalized, i.e., proportional to exp[±const.L"], with a a positive 
constant). cj 

The crucial equations for calculating the renormalized velocity are Eqs. ( |3l"| ) and (|32|). Substituting vi by zero in 
these equations and performing the Fourier transformation in Eq. ( |3l"| ) we obtain the two equations 

d m {u) = (4 - d - 2C - z) k,H - zud u K,{u) - ^Aj'(O) ^ (48) 

(27r) d |cA 2 + AM| 



A d A d \ / / / \ AM 
A t)) — 

(2nf AJ cA2[cA2 + AM 



- 2D,H - zuduDiiw) TTT^^W .Aar-Aa ■ n /.Ar ( 4 9) 



Note that Eqs. ( [48] ) and ( |49| ) are very .similar except for the trivial scaling parts. This is the consequence of the 
quantum fluctuation-dissipation theoremEJ (FDT), 



hiIm[Rf(A,oj)] ^Tiiu 
n — coth — — , 

|i?>(A,c)| 2 2rr 



ki(lo) = - 1 1 y \ n coth -±- , (50) 



which is valid in the quasi-equilibrium situation produced by the long time scales associated with the creep motion. 
If the relation (^) is satisfied in the bare system (/ = 0), then Eqs. (fl8|) and (^) guarantee its validity for any /. 
Strictly speaking the FDT is not applicable for u / 0, but its appearance here is understood to be consistent with 
our assumption of being close to local equilibrium and our use of the approximation v = in the RG equations for 
L < Lf. Using Eqs. (25), (|36|), and ( |50| ) one can see that the expressions (Eq) and (^9|) are identical up to scaling 
terms. 

Substituting Eq. (ph into (EjJ) we obtain 



ant \ ,n/l a n f \ j A * A Ad A ''(+°) D M 

= - ^AH + ^ (cA2)3 Tl [cA2 + AM] - (51) 



Note that Ti and i£; (A, t = 0) are functions of I via their dependence on Di(u). Using Eqs. (|2J), (|25j), and ( |5C| ) we 
can write T; in the form 



30 ft/Im [A(w)] coth(fyu;/2T;) 



(2n) d J 2vr | C A 2 + AM| 2 



(52) 



In a nearly static system hi and temperature Ti renormalize to zero in a trivial way, see Eqs. ( |33| ) and ([34]), and the 
elasticity c is not renormalized at all, d\c = 0. This is a consequence of the statistical tilt symmetry (it — > it + b • z, 
with b an arbitrary constant vector) of the action ([?]). 

We now concentrate on the zero temperature limit where the correlation function T; = C{" (A, t — 0) can be written 
in the form (cf., and ©)) 
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/• do; 2ft, 







This is a consequence of the identityS 

oo oo 

y dujR>(A,iuj) = J dujIm[R>(A,uj)} (54) 



following from the analytic structure of the correlator Ri as implied by causality. It is convenient to work within 
the imaginary time formalism and substitute lu — > iui. This substitution transforms the dynamic spectrum into 
Di{iuj) = rjbj + ^2 k>2 T]^ k 'u k with the advantage that Di(itv) is real and nonnegative for positive ui. The calculation 
of the velocity v of the driven elastic manifold reduces to the problem of solving the equation 

diDi(ioj) — 2Di(iw) — ziod^D^iuo) (55) 
A d A d A d A/j+of AM 
(2^) d (cA 2 ) 3 T ; [cAHAH]' 

with the initial condition D$(iuj) = r}\u>\ + puj 2 (valid for any lu) and calculating the renormalized low-frequency 
viscosity rji = d^D^iuj), lu — > 0; having found r\\ we can determine the velocity v using the relation (fi"7|). Note that if 
Dq(ilo) = Dq(—ilu) then this property is preserved under the RG transformation. 

On very short scales I < l c the disorder-dependent term on the right-hand side of (|55|) can be neglected. On 
intermediate scales l c < I < If we can substitute A/(+0) by its fixed-point value which is related to the bare potential 
through (see appendix |b|) 

A*'(+0) 2 * ee(- 2 ^^A(0) ~ ^ee^, (56) 
where we have used the relation A(0) ~ c 2 ^ 2 (A/e ic ) e in the last equation. 



C. Naive RG and "Localization" Transition 



Before embarking on the complete analysis it is instructive to see what happens if we simply keep the leading 
low-frequency form of D(iu>) as is conventionally done in dynamic renormalization group calculations. To do this 
we substitute the Ansatz Di(iu) — T)i\u)\ for all frequencies into (|5^). Calculating the integral we easily obtain 
Ti ~ {j%iA d /rji) In (^jw/cA 2 ), where we have introduced a high-frequency cutoff Q. Substituting this expression into 
( p5| ) and approximating 

AM/(cA 2 + Di(iu)) -► Di(iu)/cA 2 = muj/cA 2 , (57) 

i.e., assuming that the low-frequency asymptotics Di(iuj) — r\ito is valid for any lo and I, we obtain the equation for 
the renormalized friction coefficient 

dim-^e-^W-rii, (58) 
nrj 

with S v ~ rj^Lf.. Unfortunately, the behavior of this equation is pathological: one can see that it would imply a 
divergence of rji at a finite length scale. This would presumably mean that the velocity goes to zero in the presence of 
a nonzero force and quantum fluctuations, a result that appears to be implausible. Of course, what one must check 
in any situation where some parameter in an effective action diverges under the RG flow is, whether this is due to an 
unphysical restriction of the space of relevant parameters, a breakdown in whatever approximations that have been 
made in deriving the RG flow, or some other effect. In our case, it will turn out that the renormalization of the whole 
frequency spectrum is very important. This should be contrasted with_the classical case for which one can obtain the 
result v oc exp[— U(f)/T] by considering the low-frequency limit onlyll3. 
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D. RG Analysis of Dynamic Response 



We now turn to the analysis of the RG flow equations ( |55[ ) and (g6|). In order to find the flow of the function 
Di(uj) it is convenient to make the unconventional choice z(l) = and allow r]\ to change arbitrarily. If one is more 
comfortable with a flowing dynamical exponent z(l), one can work more generally with the quantity 

Ei(Q) = Di(m)/cA 2 , (59) 



where = wexp[— J dl'z(l')] represents the unrenormalized frequency; substituting ( p9| ) into (55) and accounting 
for (|56l ) as well as the trivial renormalization of hi , see (|34|) , we arrive at 



i -t- r. 1 1 i 1 1 i ' 



l + EASl) f +oc ndn/cA? 

J C 



2J? ' ( " )+ i + g l( n) (60) 



with the dimensionless constant K given by 



and 



K = 2Tre-^-e-( d+2 <-^ ~ §I e -(W= ( 6 1) 
(cA z ) z h n 



Ke (d+2 ( -2)l V ( L ) 



J-oo 1+B,(0) J-oo ! + £,(£!) 



For 2 < Z c the renormalization of the dynamics due to disorder can be neglected, as properly expressed by the step 
function 0(1 — l c ) in (^0|). The quantity Bi cx governs the fluctuation-induced smearing of the cusp in the 

force-force correlator. The right-hand side of ( |60| ) behaves differently for Ei(Ci) <C 1 and Ei(Cl) ^> 1; using the 
approximation 

El{n) ~min (£?!(").!) (63) 



considerably simplifies the analysis but does not change the result qualitatively. We rewrite (60) in the corresponding 
(2-dependent) frequency regions in the form 

n 2^(0) + Ei(Q)B t , Ei(Q) < 1, (64) 
9; m 2^(^) + B;, > 1. (65) 



As is readily seen from (|64|) and ( p5|) the function Ei(Q) is an increasing function of I for all fl and Ei(Q,) — > oo as 
2 — ► oo . In addition, if i?o(^) is a monotonically increasing function of Q then Ei(Q) remains monotonic in f2 for any 
I. Let us define the frequency fi; and the scale In which solve the equation 

E l (Q) = 1; (66) 

the function f2; starts at I = with a finite value fio and decreases with increasing I. For any I > we distinguish 
between the three regions < fl < f2;, f); < SI < f2g, and Jl < see Fig. ||. 
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FIG. 3. Different regimes relevant in the integration of the dynamical equation (|30|): For Q < Clo the integration is split into 
two regimes < I < In and In < I, while for > Qo the integration involves only one regime. 

In each of these regions the function Ei(n) can be found explicitly in terms of Eo(Q) and the yet-to-be-determined 
function Bi, 



Ei(n) = E (n)exp 



21 



t 

J dl'Bv 



1+ / dl'e-^'-^B,, 



;(Q) 
i 



Ei(Q) =E (tt)e M +e M J dl'e' M B v , n <Q. 
o 



, 0<n<Qi, (67) 
£li < Q < n , (68) 

(69) 



Note that we integrate (|64|) and (|65|) subsequently in the first two regions < I < Iq and In < I, while for S7 > Qq 
equation ( |65| ) applies for all values of I. Correspondingly the integral /_ ™ dSl/[l + Ei(Q)] that determines Bi, see 
(^32|), can be written as a sum of three terms. We will show below that Bi increases exponentially, implying that the 
boundary fi; is exponentially small and the first integral extending over the interval < £1 < can be neglected, 
hence 



dQ 



1 



l + Ei(Q) 



dn 



+ DC 



I + e 2(z-;(n)) 



1+ / dl'e- 2 (!'-~ l W)B v 

dn 



dn 



4 1 + E Q (n)e 21 + e 2i / dl>e-WB v 
o 



1 + G, '7 l + £ (tt)e 2; + Gf 

In the last equation we have approximated SI/ ss and Z (S^2) ~ and have introduced the expression 

l 



(70) 



(71) 



Gi = e 21 J dl'B v e 
o 



-2/' 



dhxBi/dl- 2' 



(72) 



the last approximation applies if -B; increases faster than e 21 . Substituting E (n) — (r/n + pfi 2 ) /cA 2 (see (38) and 
(|59|)) into ( |7l| ) we see that the second integral on the right-hand side of (71) always dominates and we conclude that 
it is the high-frequency behavior of Dq(uj) that controls the large length scale renormalization of the dynamics. 
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Using the definition of Bi , see 
Bi valid at large length scales 



and evaluating the integrals over f2 explicitly we obtain an implicit equation for 



Bt 



Ke (d+2Q-2)l 



„2l„2 



-21 l n 



v 2 /p e -i v/Vp 



G,cA 2 



VG^A 2 



(73) 



where we have ignored multiplicative factors of order unity. The first term in the denominator applies for small 
p at intermediate length scales when G; is not too large, while the second is relevant if the inertia is substantial 
(p > n 2 e 21 /cA 2 Gi) and for asymptotically large scales for any non-zero p. Given that the approximation in (|7^) 
is valid, as is the case in the regimes of interest, we have a non-linear differential equation for Bf y in its simplest 
approximation with G/ ~ Bi this reduces to an algebraic equation which is readily solved. Using the definition of K , 



Eq. (pl|) and the expressions L c = A 



and S 2 



p£ U C L^, we find the result 



Bt 



p cA 2 



R 2 |- e (d+2C-l)/j 2 



p_ e {d+2i-i){i-i c ) 



(74) 



in the inertial case. The behavior is rather more complicated for the dissipative case: at intermediate length scales 
we have 



B\ 



3 (d+2C)(*-U 



\v.[{r] 2 / pcA 2 K)e-( d+2C .-V l ] h ln[(S„fi/ 'Sa) e -((H-ac-2)(l-i e )] 



(75) 



which increases slightly faster with length scale than in the absence of the logarithmic factor. The crossover between 
the dissipative and massive results appears at 



■r.2 <?2 

11 pcA 2 S 2 



comparing with ([75]) this translates into the length scale 



with a corresponding energy scale 



Ui 



l/(d-2+20 



(ci+2C-2) 



hi] 



c S 2 p 



(76) 



(77) 



(78) 



For L > Lj the behavior of Bi is always dominated by the inertia and the result (Q) takes over. 

In the end we see that the coefficient Bi oc 1/r^ describing the fluctuations rounding the cusp in the correlator A; 
increases dramatically with increasing scale / (with a correspondingly decreasing Ti). Substituting E\ — * E^e 21 + G/ — > 
Bi + (Q/Q c ) a in the integral of (p2|), we see that the integration is squeezed to the high energy side where it is ultimately 
cutoff by the inertial term (a = 2) or a more general cutoff (n/fl c ) a . Hence the remaining high-frequency fluctuations 
measured with respect to the typical barriers U(L) at this scale determine the smoothing coefficient T;. Technically, 
the exponent in the non-Arrhenius type law then appears via solution of the implicit equation for Bi with the 
result Bi oc l/h a . For the extreme case with a linear spectrum sharply cut at the frequency f2 c , Dq{ilo < ilu c ) = r/uj 
and Do(iuj > iuj c ) — oo, a similar calculation provides an action that is exponentially (rather than power-law) enhanced 
in the limit of small forces /. 

Substituting the expressions for Bi into the equation for Ei we can find how the dynamic spectrum Di(iuj) is 
renormalized. At low frequencies 



and hence Di(iuj — > 0) tx exp 



diDi(iu -> iO) - [2 + B{l)]Di{iuj -> iO), (79) 
21 + j Q all' B i> ; consequently the renormalized viscosity r)f on the scale If is 

rn, ~ »yexp 21 f + / dl'B v (80) 
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(note that if we had used the conventional normalization in the action, adjusting the dynamic exponent z(l) to keep 
rji fixed, we would have obtained the same results for physical quantities but the renormalization would have gone 
into z(l) — 2 + Bi rather than the friction coefficient rj{). 



D,{m) 



cA 2 



W T1P C0 <*> 

FIG. 4. Schematic renormalization of the dynamic spectrum Di(iu>) for a dissipative/massive initial dynamics, 
Da(iuj) = t]\u}\ + puj 2 , with a crossover at uj vp — r\jp. The spectrum is renormalized differently on small and large fre- 
quencies: in the region D\(iu3) <C cA 2 the spectrum is double-exponentially enhanced and grows as t]\w\ exp[const.e' 3! ], with 
ft > an exponent depending on the specific dynamics. For Di(ico) 3> cA 2 the renormalization takes the different form 
Di(iuj) ~ tjwe' 2-2 ' 1 + puj 2 e^ 2 ~ 2z ^ 1 + Fi, with Fi independent of the frequency uj. 

Going back to the dynamical equation (|60|), we see that Bi not only determines the low- frequency part but the 
entire function Di(iuj) (we assume z = and identify SI with uj). A schematic renormalization of the dynamic 
spectrum Di(iuj) is shown in Fig. ^: the low-frequency part at u> < uoi ~ {ch 2 /if) exp(— B{) where Di < cA 2 remains 
linear, Di(iuj) ~ r/iuj, but is boosted exponentially with rji ~ rjexp(Bi) (with Bi itself growing exponentially in I). 
At frequencies above uji the upward renormalization is reduced and the response Di{iu>) develops a flat intermediate 
regime. Finally, at high frequencies uj > ujq the renormalization remains small, while the shape of the dynamical 
response again reflects the form of the original bare dynamics Do with an additional shift oc cA 2 Bi. It is this high- 
frequency part of Di that dominates the important renormalization of rji at low frequencies. We attribute the strongly 
renormalized low-frequency part D\ < cA 2 to those degrees of freedom of the manifold describing its inter-valley 
motion, while the remaining modes at intermediate and high frequencies describe its intra-valley motion smoothing 
the disorder landscape. In the next section we will discuss the meaning, significance, and problematic aspects of the 
unusual dynamical renormalization scheme uncovered above. 




CO/ 



E. Results: Quantum Creep 



The physical quantities of primary interest can be obtained from the analysis of the previous section. In particular, 
using the relation (E7h we find the creep velocity at low forces in the inertial case, 



v oc exp 



■% ffc 



(d+2C_l)/(2-C) 



(81) 



with an unknown multiplicative coefficient of order unity incorporated into S p . In the dissipative case there is a 
characteristic crossover force 



// ~ fc 



pUc 
hi] 



2/(d-2+2C) 



that separates two distinct regimes: For intermediate forces // ^ f f c we have 



v oc exp 



Sjj / fc_ 

' n V/ 



(d+2C)/(2-C) 



1 



Hf/fi) 



(82) 



(83) 



with an unknown constant coefficient incorporated into S v ■ But at asymptotically low forces / <g ; // the inertial term 
in the action dominates and the behavior of the velocity crosses over to the inertial result (pi]). For the dissipative 
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case, the result is similar to that obtained using scaling arguments in section [II B in the intermediate regime, but 
the logarithmic factor causes the velocity to decrease slightly faster with decreasing force than anticipated. By 
contrast, in the inertial case the creep velocity is much smaller that anticipated: i t can be written in the form 
v oc exp[— (S mass (f)/h) 1, with S mass (f) the characteristic action obtained in section [II C using scaling arguments. 
At asymptotically low forces, this result is also valid for the more general model including a dissipative dynamics at 
low frequencies and an inertial dynamics at high frequencies. 



F. Results: Classical Creep and Crossover 

At high temperatures the coefficient T/ = C(" (A, t = 0) is independent of the dynamics and the evolution of the 
spectrum Di(ioj) does not feed back into T;. Equation (|79|), written in terms of the bare temperature, then takes the 
form 

d l D,(iu^0)~!^D l (w^Q). (84) 

Using the relations @, U(L) ~ U c (L/L c )( d+2 ^ 2 \ L ~ L c (/// c ) 1 /(2-C) ) we obtain the classical creep law v{F) ~ 
exp[-U(f)/T] with U(f) ~ [/ c (/ c //) (d+2C ~ 2)/(2 ~ C) , see Ref. || The exponent 6 = d + 2( - 2 that determines the 
scaling of U (L) is simply the scaling dimensionality of the temperature, see (|33|), in terms of which (|84|) can be written 
with a prefactor U(L)/T ~ U c /Ti ~ e^+^-^'^Uc/T. 

One can find the crossover temperature from classical to quantum creep by comparing the exponents in ( |83| ) and 
(H|) with U(f)/T. For the inertial case, 

T cr cx ^+20/(2-0 (85) 



which is different from the naive result obtained via simple scaling arguments, see section [II. For the dissipative case, 
the crossover temperature depends on the regime. For intermediate forces, it is the same as that given by the naive 
scaling arguments (up to a logarithmic factor) , but for / <C fj the crossover temperature is again given by the result 
( |S5|) for the inertial case. 



G. Interpretation 

As mentioned earlier, a dependence on h of the form ( |si| ) might be expected if the dynamics were dominated 
by atypical barriers. A simple example is given by a classical particle diffusing in a short-range correlated random 
potential. At positive temperatures there is a linear response to an external force with the inverse mobility of the 
particle proportional to J d\JV{U) exp(U/T), where V{U)d\J is the probability density of barriers U. For the case 
of a Gaussian distribution of the random potential, we see that v oc / exp(— const. /T 2 ), resulting in a non-Arrhenius 
temperature dependence. 

Analogous effects can occur in quantum transport: Consider a quantum particle of mass m and with a dissipative 
coefficient rji, tunneling through a succession of barriers of random heights U but, for simplicity, uniform width a. 
The inverse mobility of such a particle can be written as 

- oc J dUV(U)expl-^- + — - — J. (86) 

The exponent contains a sum of two actions, the first due to dissipation and the second describing the inertial 
response. One can see that even if the dissipation is strong, in the limit h — > the inertial effects can dominate: 
the integral in (^) is calculated using the method of steepest descent and since the massive contribution to the 
action is proportional to y/U, it will contribute a term larger than the dissipative one. In particular, with a Gaussian 

distribution of barriers of the form V(U) ~ exp|— (U/Uq) 2 ] this leads to fi oc exp[— (^mUoa/h) 4 ^ 3 ] and we obtain a 
similar non-trivial dependence on Ti as found above in ( J8l| ) for the creeping elastic manifolds of interest here. 

In this simple single-particle example it is easy to understand what is going on: Because the particle must tunnel 
through a succession of barriers, the dynamics is dominated by the largest ones as long as there is a sufficiently long 
tail to the barrier distribution; the smaller the quantum fluctuations, the larger the barriers that dominate. For the 

probability distribution function V(U) chosen above, U ~ (\fmaUQ /fi) 2 ^ — > oo for ?i — > 0. The form of the tail of 
the distribution of barriers thus dominates the mobility and gives rise to the unusual dependence on h. The fact that 
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for small Ti the dominant barriers are high implies that the characteristic time scale for tunneling through the barrier 
is short (in the above case the time r is given through r ~ a^J m/U); this is what causes the long-time behavior as 
manifested in the mobility to be dominated by the inertia rather than by the dissipative response. In the case of 
interest here, the elastic manifold, the barriers that must be surmounted depend on the driving force — the lower 
the force, the higher these barriers are. If one assumes that the barriers relevant for tunneling have an unbounded 
distribution similar to the toy model studied above, it is not surprising that it is the inertial dynamics that dominates 
the limit of low forces. 

However, the reason for the anomalous dependence on h is more subtle for the elastic manifold with its many 
degrees of freedom. As can be seen from the analysis in the previous section, the anomalous dependence can be 
traced back to the dependence of the quantum fluctuations on one length scale, as parametrized by Ti, on the random 
pinning at smaller scales. Very crudely, this might be interpreted as leading to an increase of the effective mass 
density with length scale caused by the motion of smaller scale sections of the manifold implicit in the tunneling 
motion of a segment of size L f . Determining whether or not this is a reasonable interpretation must wait for a better 
understanding of the physics underlying the RG results and whether these are valid. 



V. VALIDITY OF RG RESULTS 



In both the previous and the present work on classical and quantum creep of elastic manifolds the validity of the 
approximations that underly the RG formulation have not been carefully examined. In unpublished worked one of 
the problematic aspects, the possible effects of tails in the distribution of local effective friction coefficients, has been 
investigated. Here, we briefly summarize the potential problems that this suggests as well as more basic ones that 
have not, to our knowledge, been raised previously. 



A. Random Friction 



One difficulty, analyzed in Ref. is already apparent in the toy model of a single particle in a rando m potential 



the broad distribution of times to go through or over barriers. In particular, as discussed in section IV E above, for a 
single particle the mobility is dominated not by the typical or even the average rate for overcoming the local barriers, 
but by the average time to overcome them; and the average time is dominated by anomalously large barriers. As this 
problem already arises in the classical case both in the toy model and for elastic manifolds, we restrict our discussion 
to the simpler classical limit. 

The main idea of the RG is to derive equations which relate the renormalized parameters of the field theory to the 
bare ones. Very often the parametric space of the bare and renormalized theories are identical. In other words, if the 
bare theory is described by the parameter set o.q~\ . . . the renormalized theory will be described by the same set 
of variables , af^ .... An example of this kind of RG is the </> 4 theory to one-loop order. It could be, however, that 

under the RG flow additional variables /?j , . . . are generated even if their bare values are zero. These variables 
might be strongly relevant and feed back to the original set of parameters. An example of such a behavior, which we 
could handle successfully, has already been considered above: in order to obtain sensible results we had to introduce 
a function Di(iu>) describing the dynamics on all frequencies. In this paragraph we will show that another set of 
dangerous variables is generated under the RG flow — these variables describe the probability distribution function 
of waiting times. In the RG scheme considered above the most crucial quantity is the renormalized viscosity r)i which 
is proportional to the waiting time at the scale I. We will show that the randomness due to the point-like disorder will 
produce a random and spatially inhomogeneous distribution of frictions which appears to be very broad and hence 
cannot be properly described by its first moment r\i alone. 

Let us then consider randomness in the local effective friction coefficients from the beginning and assume that the 
friction r) is a function r\ — r][u(z),z\ of both the displacement u(z) and the internal coordinate z of the manifold. 
The local rj has the natural interpretation as a characteristic time to overcome barriers involving the smaller length 
scale deformations that have already been integrated out. The simplest case to consider is a random potential that is 
periodic in u with a locally random phase shift; such a model is applicable to charge density waves (CDWs). Because 
of the periodicity, for CDWs we can expand the function rj(u, z) into the Fourier series and there will be a component 
that is independent of u and only depends on z; we consider the effects of such randomness here. 

We assume that r](z) is a random short-range correlated variable with cumulants rf n >] e.g., the first three take the 
form t](z) = = 77, r)(z)r](z') — r](z) r)(z') = ?/ 2 )<5(z — z'), and 
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ry(z)r/(z')r/(z") — r](z)r)(z') r/(z") — r)(z")r)(z') rj(z) — r)(z")rj(z) rj(z') + 2rj(z) 

= 2r/ 3) [5 (z - z') <J (z - z") + <5 (z' - z") <S (z' - z) + 5 (z" - z) 5 (z" - z')] . (87) 

Note that we define cumulants rf n > up to a factor n\. After averaging over the randomness the classical MSR-action 
(0) will have additional terms of the form 

Aand = A(n) = (-l)V n) / d d zdh . . . dt n «(*, h)iy{z, h)... u{z, t n )iy{z, t n ). 

n>2 n>2 

When deriving the RG equations with an action that includes terms of the form A ran d, it is necessary to find the 
average over fast modes of terms containing products of two perturbations, in particular, terms of the form 



5 A = i /A ra nd J d d ZTidT 2 Ai [li(z,Tl) - u(z, 7>j)] iy(z, Tl)iy(z, T2, 



(89) 

with (•••)> the standard RG average over fast modes. Making use of (|B8|), the average ( |89|) can be written as a sum of 
terms involving the cumulants 77, (see appendix^ for details; here, we summarize the main ideas of the calculation). 



The term of order n in (89) then generates 2n terms proportional to SA^\ 

8A^ ] = (-l) n vl n) A{'(0) J d d zdt l . . . dtn u{z, h)iy{z, h) . . . u(z, t n )iy(z, t n ) (90) 

which feed back to the term (88). In addition, there are An 2 — An terms proportional to 5A% , 



8A^ ] = (-l) n Vi n) J d d zdt 1 ...dt n ii{z,t 1 )iy{z,t 1 )... 

zi(z, t n )iy(z, t n )Ai [u(z, ti) - u(z, t 2 )\ ■ (91) 

These renormalizations involve the behavior of the correlator A;(u) at the origin whose growth at long length scales 
is crucial for the renormalization group analysis, see the above discussion. We can keep track of the most dangerous 
terms by substituting A" (u) by A" (0) in (j^j . The renormalization of the cumulants rfc from the above process can 
be written in the form 

d4 n) oc-A^O)^ 2 -™)^; (92) 

furthermore, there are non-linear terms that create higher moments from lower moments, see appendix As — A"(0) 
grows exponentially with Z, i.e., as a power of the length scale, the cumulants r/^ grow very rapidly. The n 2 coefficient 
and the positivity of — A"(0) oc 1/2] cx e in Eq. (92) imply that the high order moments grow so fast that ratios of 

the form m /rff, which naively are expected to be dimensionsless (in the RG sense), themselves grow exponentially 
with increasing length scale. Indeed, the high order moments increase so rapidly with n that, if these results are taken 
literally, the distribution of rji (z) has such a long tail that it is not uniquely determined by its moments — and it is 
certainly not well characterized by its mean 77; (z) = rji. Note that a random friction 77; (z) is not dangerous near the 
zero-temperature depinning transition (/ — f c <C f c ) as in this case A"(0) should be substituted by AJ'(0+) > and 
Eq. (^) suggests that -q^ renormalizes to zero for any n (although it will actually be stabilized at a small value of 
the order of an n-dependent power of e because of other terms). 

The analysis in appendix |C] shows that even if initially the friction is non-random, the disorder term alone will 
generate the corrections to the second cumulant. The second cumulant will then generate the third cumulant to 
the next loop order, and so forth. As the cumulants grow extremely rapidly, the RG flow becomes essentially 
uncontrollable. However note, that there is still an approximation in the above analysis: we have substituted the 
argument of the second derivative of the disordered correlator A"(u) by zero in Eq. (pi|). In order to be accurate we 
have to include the terms of the form 



= {-l) n r)\ n) J d d zdt 1 ...dt n ii{z,t 1 )iy{z,t 1 ). 



ii(z, t n )iy(z, t n ) Y F i [ u { z > ~ u ( z > tj)] > ( 93 ) 
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cf. Eq. (J9l|), into the action, i.e., we have to renormalize one more function Fi(u). Under the RG transformation 
this function will produce a new set functions, and so forth — it is presently unclear how all these variables can be 
analyzed in a regular way. 

One way that one might hope to make progress is to rewrite the equations of motion so that instead of having 
to deal with a random r\ on the left hand side, one works with a random mobility \i = 1/rj on the right hand side. 
This quantity, as it is bounded from above by the fastest motion, is unlikely to have troublesome long tails in its 
distribution. But the appearance of a random mobility multiplying all the spatially random pinning forces introduces 
additional technical complications into the formalism and it is presently not clear how to handle them. Nevertheless, 
there are a lot of constraints on the renormalization, e.g., the static response will not be modified by the randomness 
of either the pinning or the mobility at any wavelength. Whether this is enough to make possible a fully controlled 
analysis of at least the classical thermal creep regime is an interesting challenge. 



B. Underlying Formulation 

It is possible that the apparent runaway of the distribution of the local friction coefficients is indicative of a 
breakdown in the basic scaling assumption that underlies this and earlier work: If the barriers for motion scale with 
an exponent ip that is larger than the exponent 9 = d—2 + 2Q controlling the scaling of fluctuations in minimal energies 
and the renormalization of the inverse temperature, then the present scheme where the dynamics is controlled by static 
properties, such as the correlator A(u), cannot be valid. One then has to take into account all the dangerous variables 
discussed above. 

It is instructive to go back to the original formulation of the RG expansion for the depinning in the absence of 
fluctuations^, henceforth NF. In the derivation of the </> 4 theory from the Ising model for conventional equilibrium 
phase transitions, the starting point is an expansion around a mean field theory and the actual "field" </>(x) used in 
the RG formulation is closely related to the local effective field — applied plus exchange — acting on a spin rather 
than to the spin itself. If the interactions are long but finite range, these fields will be slowly varying in space and 
weakly fluctuating, enabling a systematic expansion to be started. 

NF focus on one segment z and use the linear combination u(z,t) of the displacements of other segments that 
determine the elastic force on z as the basic field, which is hence intrinsically a coarse-grained quantity. The underlying 
local displacements we will here call w(z, t). The segment z feels a linear restoring force proportional to u(z, t) — w(z, t), 
plus the applied driving force, plus a quenched random pinning force that is a function of w(z), and thermal noise. 
The vertices in the effective field theory are given by correlations and responses of w(z, t) to the time-dependent fields 
u(z,t). In particular, the force-force correlator A(w) that plays an essential role is related to the average 

T=([w(z,t)-w(z,t')] 2 ) (94) 

over the random pinning forces and, at positive temperature, thermal noise. In general, T is a functional of u(z, r) 
over all times r. For zero temperature depinning, the case of primary interest, the possible fields u are limited to 
those that are non-decreasing in time. In this case, it can be seen that the crucial parts of T (which are sufficient to 
analyze the depinning critical behavior) depend only on u(z,t) — u(z,t') and the functional T simplifies to a function 
of this one variable: it is then of the form assumed for A(it). In particular, if u does not change between t and t' , w 
will not change cither unless the w at the earlier time was the cause of a jump out of a formerly stable configuration 
into another; this jumping case can be handled by putting in time delays into the definition of u and, beyond this, 
the local dynamics will be independent on the history of u prior to times t and t' . 

As soon as one considers a more general dynamics — e.g., still fluctuationless but with the applied force allowed 
to decrease with time, non-monotonic stress transfer kernels, or thermal noise — the simplification of the functional 
T does not occur. In general, it is then not clear how to proceed. In the case of interest for the present paper, 
one could first assume, as in all expansions about a mean field theory, that the fields are slowly varying in space 
and time and weakly fluctuating about a uniformly advancing solution which has a slow mean velocity v. The local 
displacements will lag behind due to the pinning but will be pulled ahead by the driving force: the balance of these 
effects determines the velocity-force relation. With thermal fluctuations, the displacements will lag less than they 
would otherwise because at low velocities they have the time to surmount energy barriers. Consider now the effects of 
a time dependent change in u(z, t) on (w(z, t)). If this change is very slow, w will follow approximately adiabatically. 
But if the change is relatively fast — as it can be due to the fast motion of a neighboring segment over a barrier — 
how the local displacement w will respond will depend crucially on whether it has already been near a surmountable 
barrier for some time, and thus is likely to have already surmounted it, or has recently arrived near a barrier and 
could thus be pushed over it by the change in u. Thus we expect the responses and correlations of w to depend on 
the whole prior (and intervening) history of the fields u. 
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The basic issue that must be addressed to make further progress is whether or not the essential information about 
the basic activation processes can be subsumed in simplified functions, such as the force- force correlator A(u — u') 7 that 
appear in the formulation used in the present paper (note, however, that as we are performing an epsilon expansion, 
the central limit theorem is likely to be helpful: as we go from one scale to another, we can always assume that there 
is a large number of segments (or effective intermediate scale segments) whose dynamics contributes to the fields on 
larger scales). One possible route to proceed is to start by considering the nature of the typical fields that will arise 
from the thermally activated motion of large segments on the dynamics at smaller scales and ask whether the local 
responses to these are typically simple enough to be captured by the type of approximations used here. 

We leave further work on these points as a challenge for the future. It is worth noting, however, that the physical 
origin of the problems, the intrinsic history dependence of the phenomena of interest, arise in many other non- 
equilibrium systems with many degrees of freedom. Better understanding of them in this context is thus likely to 
bear fruit more generally. 

C. Alternative Expansions 

Many of the difficulties encountered in this paper would be lessened if the energy scale — and hence the barriers 
— did not grow rapidly with length scale. One way to get around this difficulty might be to consider an expansion 
about a different limit, one in which not only the roughness exponent C is small, but also the scaling exponent 9 for 
the energy. In d-dimensional manifolds with long range elastic interactions that fall off as i/ r d+a with < a < 2, the 
upper critical dimension for both fluctuationless depinning and rough manifolds in equilibrium is d c — 2a (i.e., £ = 
at d = d c ). The energy scaling exponent is 9 — d — a + 2£ so that for a small, 9(d c ) = a is also small. One might then 
hope to attempt an expansion in a and e = d c — d simultaneously. This would have the advantage that, assuming the 
barrier scaling exponent ip is also small for small a near d c , the exponential dependence of the velocity on the length 
dependent barriers would be relatively weak and hence, perhaps, systematically controllable. We must also leave this 
and other possible limits about which more controllable expansions might be attempted for future research. 

D. Quasiclassical Langevin Equation 

In addition to potential problems with the renormalization group formulation that are analogous to those discussed 
above for classical creep, in order to obtain results for quantum creep we have resorted to a quantum Langevin 
equation to approximate the quantum dynamics. We would like to be able to investigate this formalism more deeply 
and understand its limitations and regimes of validity. It is instructive to first consider simpler problems. 
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FIG. 5. Tunneling of an overdamped particle out of a ID potential well in the presence of an external force /. The quantum 
mechanical answer for the tunneling action is of the form S oc l// 2 , / — > 0. The action obtained via the QLE is determined 
only by the region xo < x < xi. Since xi — xo <C X2 — xi, the QLE overestimates tunneling effects (i.e. gives smaller action). 

The QLE can be used to analyze tunneling of a single particle if the barrier separating the metastable state from the 
lower-energy stable one has a single characteristic height and width. In this case the tunneling action obtained from 
the QLE differs from the true one only by a factor of order unity. In principle, the QLE can lead to an overestimation 
of the tunneling rate. For example, consider an overdamped particle (with viscosity r\ ) in a ID potential well U(x) 
of depth Uq with a minimum at x = and U(x) rapidly decaying for x — > ±oo. In the presence of an external force 
/, see Fig. [BJ the well becomes metastable and the particle will leave. The quantum action S can be estimated as 

S ^ 77(^2 - xi? « v{Uo/f) 2 . (95) 

This action is dominated by a region of width ~ 1/f where the potential U(x) is negligible but fx is less than the 
depth of the well: this yields S oc 1/f 2 as / — > 0. If we use the QLE approximation, the tunneling action would be 
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much smaller than the true WKB value in the limit / — > 0: since the QLE is a local approximation, the only relevant 
part of the potential would be that for x < x\, with x\ the point where U(x) — fx is maximal. In the QLE, once 
the random force that obeys the quantum fluctuation-dissipation theorem brings the particle to the point xi, the 
particle can leave the well and the region x > x\ that dominated the WKB action does not matter. Since in the limit 
/ — ► 0, X\ — xq -C x<2 — x%, the tunneling action obtained within such an approach is much smaller than the correct 
one, Sqle — vi x i ~ x o) 2 ^ S\ the QLE can thus drastically overestimate tunneling rates. 

We believe, however, that for the quantum creep problem considered in this paper the QLE should be a reasonably 
accurate approximation for the action whose exponential determines the velocity (|8l|): in the naive scaling arguments, 
the effective potential wells that segments must tunnel out of have only one characteristic width and height. Neverthe- 
less, the results we obtain from the QLE and RG analysis yield a much larger effective action for the tunneling than 
the naive scaling arguments. The source of this is unlikely to be due to the QLE approximation as, by analogy with 
the single particle case above, we would expect that, if anything, then the QLE would underestimate the tunneling 
action. This argument is reinforced by noting that the scaling analysis in section III can be applied directly to a 
system driven by a stochastic force obeying the quantum FDT for which the QLE is exact: these would lead to the 
same results as those obtained from the scaling arguments for the full quantum dynamical problem, see also Ref. |2(1 



VI. CONCLUSIONS 



In this paper we have considered the motion of an elastic manifold driven through a disordered medium at low 
temperatures where quantum fluctuations are important. We have focused on the limit of small driving forces (/ <C f c ) 
for which the average velocity of the manifold is small and dominated by quantum tunneling through barriers between 
locally stable configurations. Using an at-least-partially controlled RG expansion, we find that the resulting creep 
velocity is exponentially small in a power of 1//. While for strong dissipation we find an intermediate range of forces 
where the creep law agrees with the result of simple scaling estimates (up to logarithmic corrections, cf. ( |l7j ) and 
(|S3|)), we find that the results of the RG analysis are more complex. In particular, at asymptotically small forces 
the creep velocity has the non-trivial form ( |gl| ) , combining an underlying massive (or other high-frequency) dynamics 
with an unexpected non-Arrhenius-like dependence on Ti. The structure of the renormalization group flow is rather 
subtle and very different from that for classical creep. In particular, because of the importance of both rapid motion 
during a tunneling event and the slow overall motion, it is necessary to considep,the dynamics of the manifold at all 
frequencies u>. This is in contrast to the classical case for which Chauve et a/J13 considered only the low- frequency 
limit of the dynamics. In the quantum case such a low- frequency approximation would lead to a spurious "localization 
transition" with the average velocity of the manifold dropping to zero at a small but finite value of the driving force. 
Although future improvements based on the analysis of the complete quantum mechanical action rather than the 
quantum Langevin approximation might change some of the results presented here, we believe that the main feature, 
the importance of the whole spectrum of frequencies for the tunneling dynamics, will persist. 

In spite of the subtleties that appeared in the analysis we have carried out there are further difficulties associated 
with both the broad distribution of barriers and the underlying field theoretic formulation of expansions about mean 
field theory. Although we have not resolved all the difficulties, and thus are not sure whether the present results are 
truly systematic, we have noted the physical and formal sources of the problems in ways that we hope will help direct 
future progress. 
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APPENDIX A: DERIVATION OF THE REAL-TIME QUANTUM ACTION 

In this appendix we discuss the derivation of the action (^): Consider the classical equation of motion (|j) with 
the (5-correlated noise /th(z, i); we search for the quantum Hamiltonian that reproduces the same equation of motion 
in the classical limit. One class of possible Hamiltonians has been introduced by Caldeira and LeggettH where the 
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system under consideration (an elastic manifold in our case) is assumed to be coupled to a bath consisting of harmonic 
oscillators. The Hamiltonian of the bath can be written in the form 



bath 



= / aj; 



if(z) Mj (z) { 



2Mj(z) 



fij(z) X,-(z)- 



Ci(z) 



M»ft 2 (z) 



u(z) 



2-, 



(Al) 



where each segment of the manifold z is coupled to the system of harmonic oscillators. This Hamiltonian contributes 
to the quantum-mechanical action with 



S bath [X] + S int [ u ,X] = J dtj d d zY^ 



M,(z) 



O (z) A^-(z)- 



2 */« 
c j( z ) 



u(z) 



2i 



(A2) 



where X is the vector with components Xj describing the oscillators of the bath. The action corresponding to the 
elastic manifold Sq[u] has the form 



S [u] 



dt / d a z 



c ( du 



2\di 



p f du 
2\~di 



+ U(u, z) - fu 



(A3) 



After substituting the action S — <Sbath + £Wt + So into Eq. (|8|) we can eliminate the bath degrees of freedom appearing 
only quadratically in the action. We define the spectral density 



c?(z) 



2 



cj > 0; 



(A4) 



the ohmic kernel J(w) = r/uj produces the action (^). In the classical limit (high temperatures) one can rigorously 
expand the potential energy terms in (|9|) in y and obtain the action (^) which is equivalent to Eq. (Q). Note that we 
have introduced new coordinates u and y in the action (||) corresponding to the "center of mass" and relative motion 
of the trajectories in rf§). 



APPENDIX B: FIXED-POINT FUNCTION OF THE CORRELATOR 

In this appendix we derive the relation (|5^) . We investigate the static RG fixed-point in the absence of a driving 
force; the temperature and quantum fluctuations renormalize to zero, C, > (A, t — 0) — > 0, and hence we need to 
consider the function A;(u) alone. For I < l c we can neglect the nonlinear terms on the right-hand side of (pq). 
Solving the remaining linear equation we obtain 

A lc ~ A(0)e (e ~ 2C) ^ (Bl) 

for the correlator height and ^* w ^ exp(— ^Z c ) for its width (the latter result follows from integration of the second 
term in (p8|)). On the other hand, for I ^> l c the correlator is close to its fixed-point value. Using equation ( p8| ) for 
d[Ai(u) = and substituting u = we obtain 

(e - 2C) A /c (0) ~ (e - 2Q A*(0) ~ 7A*' 2 (+0) (B2) 

and replacing I = AdA d /c 2 A 4 and (~ e we obtain (|56|). Combining these results with |A*'(0+)| - A*(0)/£* we find 
the estimates A*(0) ~ (cA 2 ) 2 ^ 2 A- d e - 2 ^= and |A*'(0+)| - (cA 2 ) 2 £A- d e -<^. 



APPENDIX C: DANGEROUS RELEVANT VARIABLES 



As mentioned in section, y A , the RG flow generates dangerous relevant variables corresponding to the cumulants 
of the friction distributioncil. From a physical point of view, the friction coefficient "q in the equation of motion (^) is 
always spatially inhomogeneous due to the presence of disorder, so that the bare values of the cumulants describing 
the friction distribution are always nonzero. It turns out that under the RG transformation these cumulants grow 
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extremely rapidly and the description of the t herma lly activated or quantum motion in terms of the renormalized 
viscosity r\i or the function Di(u) (see section [V D| ) alone is not appropriate. Since the renormalized viscosity on 



a certain lengthscale corresponds to the waiting time on this scale, we conclude that the broad distribution of the 
friction (it is broad as the cumulants grow very rapidly) implies a broad distribution of waiting times. Hence, the 
proper description of the problem requires renormalization of the waiting time probability distribution function rather 
than its first moment alone. The main goal of this appendix is to show how the dangerous variables are renormalized 
under the RG flow. We will also show that even if initially all of them are zero, they will be generated by the random 
pinning. We will restrict ourselves to the case of high temperatures; the same variables will be generated in the 
quantum case. 

Let us assume that the friction -q in the equation of motion (Q) is a spatially inhomogeneous function of the 
coordinate; in general, 77 = rj(u,z). For simplicity, we neglect the dependence on u (this does not change the result 
qualitatively). A physical realization where this approximation is valid is the pinned charge-density- wave: the friction 
rj(u, z) is a periodic function of u and, hence, can be expanded in a Fourier series; t](z) then is the w-independent (the 
zcroth) harmonic. Below we will assume that the ry-disorder and the A-disorder are uncorrelated and short-ranged; 
in reality, these two types of disorder should be correlated. If we take these correlations into account, the dangerous 
variables will still exist and the RG flow will be similar. In this section we assume that the response function has the 
form 1/ (ck 2 — irjicu) , i.e., we do not c onsi der th e full dynamic response associated with Di(u>); this approach is valid 



in the classical case, see sections IV C and IV F . 

Averaging over the disorder in 77 within the MSR-functional (0) we see that the friction term will be transformed 
into 



d d zdt riw(z, t)iy(z, t) — > — / d d zdt rju(z, t)iy(z, t) + A ran d 

f d d zdt riu(z,t)iy(z,t) +J2 A(k) ( C1 ) 

J k>2 

d d zdt rju(z,t)iy(z,t) + ^(-l)V fe) / d d z dUii(z, U)iy(z, U). 

J ,; — 1 



k>2 

When deriving the RG equations we consider the term A ran( j and the usual disorder-induced term 

A dis = ^J d d zd Tl dT 2 A [u(z, n) - u(z, t 2 )} (C2) 

as perturbations. To one- loop order we have to calculate the second cumulant of the perturbation. There will be a 
cross-term of the form S(A Tan< iA < iis}>, where S(H} > denotes the change in H after averaging over fast modes. Let 
us show that the average of this form gives rise to the singular renormalization of the coefficients r/ n '. Figure || 
summarizes the various symbols appearing in the diagrammatic expansion: Lines without arrows denote displacement 
fields u, lines with arrows at the end correspond to the auxiliary fields iy, and lines with arrows in the middle denote 
response functions. Crosses stand for time derivatives and the wavy horizontal lines denote disorder correlators 
A; [u(z,ti) — u(z, ti)]. Dashed horizontal lines stand for (— l) k ri^ and dots . . . represent further omitted lines. For 
each response function with an arrow coming into a wavy vertex (see e.g., the diagrams in Figs. ||(a) and (b)), the 
disorder correlator A[w(z, t{) — u(z, t 2 )] is differentiated with respect to its argument. These derivatives appear as a 
consequence of averaging of the correlator and the iy- fields, see Eq. (C4) below. Dashed and wavy lines in diagrams 
connect times, the upper line connects t\, t 2 , ■ ■ ., and the lower line connects the times T\ and t%. Note that there is 
no time ordering corresponding to vertical or horizontal directions. The fields at the upper vertices in diagrams are 
taken at the spacial coordinate z and the lower vertex is at z'. The argument of the response function involves the 
differences between the times and spatial coordinates of the end and starting points. We quote three useful formulae 

5(u(z,t)iy(z',t')) > = Rf-(z-z',t-t')dl (C3) 
6 (A [u(z, n ) - u(z, t 2 )] iy(z', t')) > = A' [u(z, n ) - u(z, t 2 )] [R> (z - z', n - t') - R>(z - z', r 2 - t')]dl (C4) 

f d d z'Rf{z',t x )R?{z',t 2 ) = ^ 3 dlRi(A,t 1 )Ri(A,t 2 ), (C5) 
where i?;(A,i) denotes the (partly Fourier transformed) response function (p^) and Rf(A,t) is defined in (|3G|). 
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(a) 



(b) 



0\ 



(c) 



R 



(d) 



A 

(e) 



X 



(f) 



FIG. 6. Elements used in diagrams: (a) the displacement field it; (b) the iy-field; (c) the response function R >m , (d) crosses 
denote derivatives, with respect to time in the u-field if the cross is on a solid line, or with respect to the time difference 
in the response function if the cross is on a response function line; (e) diagrammatic representation of the correlator A; (f) 
diagrammatic representation of random friction terms (— 1)™?/™' Y\ii(z,tk)iy{z,tk), see (CI); dots . . . stand for further omitted 
lines. 



The contribution SDi of the diagram in Fig. 0(a) has the form, 



(n) 



d d zd d 



d,T\ dT2 



.k=i 



S (u(z,ti)iy(z',Tx)) > S (u(z,t 2 )iy(z',T 2 )) : 



xiy(z, h)iy(z, t 2 )A l [u (z', n) - u (z', t 2 )] u(z, t k )iy(z, t k ) 



(C6) 



fe=3 



2* 



(27T) 



d d 2 



jj dt k 



Ri{A,tx - Ti)Ri(A,t 2 ~ T 2 )iy(z,t 1 )iy(z,t 2 ) 



/n 
dridr 2 Ai [u (z,ti) - u (z,t 2 )] ]^[ u(z,t k )iy(z,t k ) 

k=3 



In the last equation, we have used (C3) in order to express averages over fast modes through response functions and 
have integrated over z' using ( |C5|) . We rewrite the product of the two response functions as i?/(A, t\— ri)i?;(A, t 2 ~ t 2 ) = 
(9 Tl i?;(A, t\ — ri)9 T2 i?;(A, t 2 — t 2 ) and integrate by parts over variables T\ and t 2 , thus generating the second derivative 
of the correlator A. The response functions connect times t\ and t\ and t 2 and t 2 . Consequently, we can set t\ ~ t\ 
and t 2 w t 2 in the time arguments of the iy-fields and integrate over the variables t\ and t 2 in the response functions 
using the formula J dtR(A,t) = R(A,oj = 0) = 1/cA 2 . In the end we obtain an expression involving the integrals 
J d d z dT\ dr 2 dt% . . . dt n and performing substitution t\ — * t\ and t 2 — > t 2 we arrive at the final contribution of the 
diagram 0(a) 



dl 



1 („) A d A d 

2 (27^ (cA 2 )^ 



d d z 



Y[ dt k ii(z,t k )iy(z,t k ) 



Lk=l 



AJ' [u (z,t 1 )-w(z,t 2 )] 



(C7) 



Note that in the diagrammatic language integrating by parts corresponds to drawing u-field lines out of the wavy 
vertex and moving crosses (i.e., time-derivatives) from the response functions down to the displacement fields. The 
low-frequency expansion of the diagram in Fig. 0(a) then takes the form of the diagram in Fig. |^(f ) [this low-frequency 
analysis involves the expansion of the argument u(z,t\) — u{z,t 2 ) of the correlator A in a Taylor series in powers 
of (ti — t 2 ); the first term (equal to zero) of this expansion produces the desired contribution]. There are n(n — 1) 
equivalent diagrams of the type in Fig. 0(a) allowed by the permutation symmetry. 



2G 




(a) 



(b) 




(c) 













( ••• : 


/ v 



(d) 

FIG. 7. Dangerous diagrams contributing to the renormalization of the cumulants of the friction-coefficient distribution. 
These terms arise from averaging over fast modes in (A ran d J 4dis}>: (a) with multiplicative factors n(n — 1) (due to the 
permutation symmetry), (b) with multiplicative factor 2n, (c) with multiplicative factor 2n(n — 1), (d) with multiplicative 
factor n(n — l)/2. For convenience we draw some of the iy-fields and response functions at a 45 degrees angle or deform the 
lines. 



The expression corresponding to the diagram in Fig. ^|(b) is given by 



SD 2 = ~p\ n) J d d zd d z' J d Tl dT 2 



Y[dt k 

k=l 



6(u(z, h)iy(z', ti))><5(A; [u (z', n) - u (z', t 2 )] iy(z, ti))> 



1 Ww„ , A d A d dl 

= — ~Tji A; (0) — -7 

2 " ' V ' (2n) d (cA 2 ) 2 



Z/( z '> T 2) \\ u(z, t k )iy(z, t k ) 

n 

d d z JJ dt k u(z, t k )iy(z, t k ); 



(C8) 



fc=2 



fe=l 



there are 2n topologically equivalent diagrams of this class. Note that in contrast to the term 8D\ the argument of 
the function A; (it) in the above expression is zero. The diagram shown in Fig. Q(c) gives the contribution 



$D 3 = ijjW J d d zd d z' J d Tl dT 2 



Y[dt k 

,fc=i 



5(u(z, h)iy(z', t 1 )) > 8(A [u (z', n) - w (z', t 2 )] iy(z, t 2 )>; 



iy(z, ti)u(z, t 2 )iy(z', r 2 ) it(z, t k )iy(z, t k ) 



(C9) 



1 (rp ^A* g 
~2 r/i (27r) d (cA 2 ) 2 



I [ dt k u(z,t k )iy(z,t k ) 



Lk=l 



AJ' [u (z,t 1 )- U (z,t 2 )]; 
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there are 2n(n — 1) topologically equivalent diagrams of this class. Finally, the diagram in Fig. |?](d) is given by the 
expression 



SL> 4 = 



(n) 



d d zd d z' I dndn 



Y[dt k 

.Ai=l 



u(z, h)u(z, t 2 )iy{z\ n)iy(z', t 2 ) 



Y\_ u(z,t k )iy(z,t k ) 



.k=3 



I (n) A d A d dl 
2^ (27r) d (cA 2 ) 2 



S(iy(z,ti)S (A[u (z',ri) - u(z' ,T 2 )]iy(z,t 2 )) > ). 



Y[dt k 



Lk=l 



dTidT 2 u(z, h)u(z, t 2 )iy(z, Tx)iy(z, t 2 ) 



! J u{z,t k )iy{z,t k ) 



.k=3 



x {R(A,Ti-ti)R(A,Ti-t 7 )+R(A,T 9 -ti)R(A,T 2 -t 3 ) 

-R(A, n - h)R(A, t 2 - t 2 ) - R(A, n - t 2 )R(A, r 2 - ti)}A" [u (z, n )-u (z, r 2 )] = 25D 3 . 



(CIO) 



The first two terms in the curly brackets do not feed back to the random friction cumulants, while the third and fourth 
terms give the desired contribution. The multiplication factor of the diagram in Fig. 0(d) is equal to n(n — l)/2. 

In order to find the contribution to 77, we need to set the argument u(z,ti) — u(z,t 2 ) in the correlator A(u) to 
zero in (C7), (|C9|) , a nd ( plQ|) (as onl y thi s term feeds back to the random friction cumulants) and sum over the four 
contributions (|C7|) , (08), (p9|), and (C1C) each multiplied with its appropriate weight; the result takes the form 



Sr)™ = -rjl n} (2n 2 - n) 



A d A d Aj'(O) 
(27r) d (cA 2 ) 2 



dl. 



(Cll) 




(b) 

FIG. 8. Diagrams describing the generation of the second cumulant of the friction distribution; their low- frequency expansions 
feed back to the second cumulant. 

Next, we show that the term A^ itself generates the second cumulant rf^ . The dangerous relevant diagrams arise 
from the average (j4di s ^4dis)>, see Fig. |[ their low-frequency expansion gives the contribution to the random friction. 
The contribution of the diagram |8|(a) takes the form 



SD i = 77 



iM d 



dl J d d zdt\dt 2 A; [u(z, tx) - u(z,t 2 )]iy(z,t 2 )iy(z, r 2 ) 
x J dridr 2 A J [u(z, n) - u(z, r 2 )] R t (A, t 2 - Tx)Ri(A, r 2 - tj) 



(C12) 



where the averaging over fast modes and the integration over z' has already been carried out. The response functions 
connect the times t 2 , n and r 2 , t\ and we can expand 



u(z, t x ) - u(z, t 2 ) = u{z, t 2 + (tx - r 2 )) - u(z, t 2 ) 

PS U(Z,T 2 ) - u(z,t 2 ) +u(z,T 2 )(tx - T 2 )\ 



(C13) 



similarly, u(z, tx) ~ u(z,t 2 ) w u(z,t 2 ) — u(z, r 2 ) + ii(z, t 2 )(rx — t 2 ). Substituting these expansions into (C12) and 
expanding the correlators A'(u) we arrive at the expression (note that / dttR(A,t) = ^/(cA 2 ) ) 



28 



SD 1 



l A d A d 
8(27r) d (cA 2 ) 4 



-dl 



d d zdt 1 dt 2 A" 2 [u(z,tx) 



u(z, t 2 )\ u{z, h)iy(z, ti)w(z, t 2 )iy(z, t 2 ); 



(C14) 



there are 4 topologically equivalent diagrams of this class. 
Diagram (b) in Fig. || gives the contribution 



SD 2 



lA d A d 



(2n)° 



-ydl I d zdt\dt 2 A" [u(z, ti) — u(z, ti)] iy{z, t\)iy(z, t 2 ) 



x J dTidT 2 Ai [u(z,n)-u(z,T 2 )]Ri(A,ri-t 1 )R l (A,T 2 -t2), 



(C15) 



where the averaging over fast modes and the integration over z' again has been performed already. The response 
functions connect the points T\ and t\ and t 2 and t 2 . Expanding the displacement fields in t\ — t± and t 2 — t 2 , 
substituting them into ( |C15|) , and expanding the correlator A(u) we see that 5D 2 = 8Di. There are 4 topologically 
equivalent diagrams of class (b). Summing up the contributions 8D\ and 8D 2 and multiplying the sum by 4 we obtain 
the contribution Sfji ' of the disorder term to the second cumulant, 



jt~( 2 ) 



A,A d 



(27T) a (OA*) 



(C16) 



Finally, we can w rite th e on e- loop RG eq uatio n for the cumulants of the random friction distribution; accounting 
for the terms Q, Q, (jc|), ( |C1C| ), and jC16|) we find the flow equations 



« ( 2 ) 



(4 - d - 2z) ril 2) + I 



[d + 2n — dn — zn) ?7 ; — (2 



(cA 2 ) 

{n) 



_^i_ A r(o),p\ 

(cA 2 ) 2 ' 
7 (cA 2 ) 4 ' V : 



(n) 



n ^ 2. 



(C17) 
(C18) 



The system of RG equations ( C17 ) and (C18) has been derived to lowest order in 4 — e. To next order, the second 
cumulant will generate the third cumulant. To third order, the third cumulant will generate the fourth cumulant, 
etc., i.e., all cumulants will be generated in the RG flow even if all of them (except for the first one) are equal to zero 
initially. From a physical point of view, the friction will always be random as the point-like impurities suppress the 
order parameter randomly and, hence, the dynamic characteristics of the medium (e.g., a superconductor) is random 
as well. Since — A"(0) oc 1/7} cx e el and the correction due to disorder grows as n 2 for large n, the random friction 
probability distribution becomes very broad and one needs to take into account all its moments and not just the 
friction r\i only. In order to obtain Eqs. (C17) and ( C18| ) we have made the approximation u(z,t\) — u(z,t 2 ) = in 
the final expressions for SDi, SD^, 8D\, and SD 2 ; in general, however, one needs to renormalize the full functional 
(|93|). During the RG procedure, other terms will be generated and it is unclear at this stage how to take all these 
terms into account in a controllable way. 
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